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Abstract 

The quantum hadro dynamics (QHD) model with minimal nucleon-meson couplings is gener- 
alized by introducing couplings of mesons to derivatives of the nucleon field in the Lagrangian 
density. This approach allows an effective description of a state-dependent in-medium interac- 
tion in the mean-field approximation. Various parametrizations for the generalized couplings are 
developed and applied to infinite nuclear matter. In this approach, scalar and vector nucleon self- 
energies depend on both density and momentum similarly as in the Dirac-Brueckner theory. The 
Schrodinger-equivalent optical potential is much less repulsive at high nucleon energies as compared 
to standard relativistic mean field models and thus agrees better with experimental findings. The 
derivative couplings in the extended model have significant effects on the properties of symmetric 
nuclear matter and neutron matter. 

PACS numbers: 21.65.+f, 21.30.Fe 
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I. INTRODUCTION 



Properties of nuclear matter and finite nuclei have been described with great success by 
relativistic mean field (RMF) models in recent years Jl], [2|, |3|, |], |5|, ||, [7], ||, These models 
provide a novel saturation mechanism, an explanation of the strong spin-orbit interaction, 
and a natural energy dependence of the optical potential. They are often referred to as 
quantum hadrodynamics (QHD) since nuclear systems are described in these quantum field 
theories in terms of interacting hadrons, i.e. nucleon and meson fields. Lagrangians have 
been proposed in simple versions about 25 years ago |IJ, and since then there have been 
many different treatments, extensions and applications. 

There are basically two approaches of QHD, which one may call microscopic and phe- 
nomenological, respectively. The microscopic method is the Dirac-Brueckner (DB) approach 



Ig, 11, 12, |m p| p| ffl O, 0, |IJ, |g, 2|, ||, |3| which tries to connect nucleon-nucleon 



scattering with finite nuclear systems via many-body theory. The effective Lagrangian has 
a simple form where nucleons and mesons couple in a minimal way. The coupling constants 
are fitted to NN scattering observables via the relativistic T-matrix approach. Then nuclear 
matter is treated in the Brueckner approximation and an effective in-medium interaction is 
extracted. The relativistic Brueckner approach successfully achieves saturation for nuclear 
matter which was not possible in non-relativistic methods with only two-body forces. Nu- 
cleon self energies Ej (i = scalar, vector, tensor, . . . ) in the DB approach are both density 
and momentum dependent. Below the Fermi momentum kp the momentum dependence is 
found to be not very strong but above kp it corresponds to the momentum dependence in 
the empirical optical potential. While the DB method is very successful in describing infinite 
nuclear matter, it is very difficult to apply to finite nuclei without drastic approximations, 
e.g. the local density approximation (LDA). In order to describe detailed properties of finite 
nuclei some adjustements have to be introduced |L4|, |4], |2f|. 

In the phenomenological approach QHD is applied directly to nuclear systems without 
reference to NN scattering. Simple many-body approximations, usually the mean field ap- 
proximation without exchange and correlation contributions, are employed. The coupling 
strengths of the mesons are directly adjusted to saturation properties of nuclear matter. In 
the original approach proposed by Walecka [l] the self energies are proportional to the corre- 
sponding densities. It was found that in this approach with linear meson-nucleon couplings 
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properties of nuclear matter are not satisfactory away from saturation (too large incom- 
pressibility and too small effective mass) and finite nuclei are not described well. In order 
to improve the description the QHD Lagrangian has to be extended to include medium and 
state dependent effects in an effective way. This has been achieved by adding non-linear and 
higher order interaction terms to the Lagrangian or assuming density- dependent couplings 
of the nucleon and meson fields. 

Non-linear models are the most prominent approach, in which cubic and quartic self- 
interactions of the cx-meson are introduced. This extension was proposed by Boguta and 



Bodmer |27| in order to lower the incompressibility of nuclear matter. Later the approach 



was extended to other meson fields ||28|| . Various parametrizations have been introduced 
which provide a good description of a whole range of nuclei |j, [F], |29], [30], [U|. Replacing 
the minimal coupling of the a-meson to the nucleon by a derivative coupling, Zimanyi and 
Moszkowski (ZM) obtained a particular non-linear cr-coupling after rescaling of the nucleon 



field ||32|| . Although this approach leads to a small incompressibility of nuclear matter, the 
effective nucleon mass turns out to be rather large, and, correspondingly, the spin-orbit 
interaction is too small. The ZM model with further extensions was studied in detail in 



33], p5| . Recently, this class of approaches has been put on a more systematic basis 
in the effective field theories (EFT) where the possible terms and their importance are 
systematically categorized [^, |37|, |359, B^, f20| . 

An alternative approach assumes a density dependence of the coupling vertices PTl 



43] . This extension is suggested by DB theory. However, in a thermodynamically consistent 
and Lorentz covariant model, the couplings cannot depend parametrically on the density. 
They have to be lorentz-scalar functionals of the field operators. Then the Euler-Lagrange 
variational equations lead to so-called rearragement contributions in the self-energies. The 



first serious parametrization for finite nuclei was developed in [[43]] where a dependence of 
the couplings on the vector density (VDD) was found to work best. It was shown that 
good results can be obtained in the whole chart of nuclei, comparable to the best non-linear 
models. Recently, a new parametrization within the VDD approach has been developed |44 



These models can be considered in the general framework of density-functional theory which 
is given here in terms of the density dependent vertex functions. This connection is rather 
close, since the density dependence is motivated by DB theory which includes exchange and 
correlation terms. 
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Generally it seems desirable to describe the structure of nuclei and nuclear reactions 
in a common framework. This implies that one also needs the momentum dependence of 
the self energies at least for momenta above the Fermi surface. Empirically these are well 
studied in Dirac phenomenology by fitting scalar and vector potentials to elastic proton 
nucleus scattering data jI5|, It is found that the optical potential becomes repulsive 
with increasing energy but levels off at about 1 GeV. Optical potentials extracted from 
standard RMF models are much too large and increase linearly at high energies. Recently, 
a reasonable description of proton-nucleus scattering and an improved optical potential has 
been obtained by multiplying the self-energies from RMF calculations of finite nuclei with 



functions that depend explicitly on the proton energy [47]. However, the consistency of the 
theory is lost and the change of the energy dependence at different medium densities is not 
properly taken into account in this approach. The momentum dependence of the interaction 
also becomes important in the description of heavy-ion collisions. Several parametrizations 
for a parametric momentum dependence of the potentials have been developed in different 



approaches g|, gjj |g |I) ||, |3|. 

So far, extensions of the RMF models concentrated on the effective density dependence 
of the in-medium interaction since they were mostly applied to the description of finite 
nuclei. If one goes beyond the mean field approximation momentum dependent self en- 
ergies will emerge naturally. Recently, it was shown, that an energy-dependence of the 
self-energies can improve the level density close to the Fermi surface of finite nuclei [S3]. In 
this RMF model, dynamical effects from the coupling of the single particle motion to col- 
lective surface vibrations are taken into account in a phenomenological approach. However, 
the energy-dependence was introduced only parametrically into the self-energies destroying 
the covariance of the theory. In order to retain the simple structure of the equations known 
from the RMF approximation one can imagine that it is desirable to modify the QHD La- 
grangian in such a way that it leads to momentum dependent self energies already in the 
mean field approximation. In this work we propose such a class of Lagrangians with gener- 
alized derivative meson-nucleon couplings. The various types of possible couplings will have 
characteristic consequences in the field equations. Here we limit ourselves to the study of 
infinite nuclear matter within the derivative coupling (DC) model in order to demonstrate 
the main effects and the limitations of the approach. The results can serve as a foundation 
for a later application of the model to finite nuclei. Although derivative couplings have been 
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investigated in simple approaches before, they were essentially cast into density-dependent 
forms, e.g. by rescaling the fields as in the ZM model. In EFT models contributions contain- 
ing derivatives of the density are included, but generally they do not lead to a momentum 
dependence of the self energies and they do not affect properties of nuclear matter. 

The paper is organized as follows. In Sec. [TJ the the Lagrangian of the DC model is 
introduced and the field equation for mesons and nucleons are derived in the mean field 
approximation. In the following section the DC model is applied to infinite nuclear matter 
and expressions for basic quantities are derived. Three parametrizations are developed 
in Sec. [IV] and the resulting properties of symmetric and asymmetric nuclear matter are 
compared to standard RMF models. Finally, we summarize and give an outlook in Sec.[V|. 



II. LAGRANGIAN DENSITY AND FIELD EQUATIONS 

The description of nuclear systems in QHD starts from a Lagrangian density £ that 
contains nucleons and mesons as degrees of freedom. In case of symmetric nuclear matter 
it is sufficient to include isoscalar a- and tu-mesons. However, for asymmetric matter and 
finite nuclei, isovector mesons have to be added. It is customary to consider only the vector 
p-meson but for the sake of completeness we also take the scalar 5-meson into account. 
The standard approach is to couple the nucleons minimally to the mesons. In case of the 
isoscalar mesons this leads to vertices of the form x/xpa and tpj^ipuj^ corresponding to the 
linear Walecka model. 

The idea of the DC model is to introduce explicit couplings of the mesons to the first 
derivative of the nucleon field if), e.g. ipid^ip or ip'y fl id fl ip. The derivate id^ in the new 
coupling vertices is actually replaced by the covariant derivative iD^ as it appears in the 
usual kinetic contribution to the Lagrangian. In the current model we consider couplings 
to the covariant derivative of the nucleon field which are linear or quadratic in the a- and 
cu-meson fields. For simplicity for the 5- and p-mesons only linear couplings are taken into 
account. Of course, the standard minimal couplings are retained in the DC model. 

On these conditions the Lagrangian density in the DC model assumes the form 



#>Z^ + (iD^)T^ -i,(m-Y a a- T s t -Pj^ + C m (1) 
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with the covariant derivative 

iD^ = idf, - TuUJft - T p t ■ p M . (2) 

The minimal coupling of the nucleons with mass m to the mesons i — a, u, 5, p of mass rrii 
is specified by coupling strengths IV The contribution 

C m = l - [d^ad^a - mla 2 + d»5 ■ dj- m 2 5 6 ■ 6 (3) 

-\g^G, v + mlu^ - ■ H, u + my ■ p, 

to the total Lagrangian density C with field tensors 

G^ u = d^u u - d u u^ = d^p v - d„Pn (4) 

describes the free meson fields. The Lagrangian density (^]) is symmetrized in the covariant 
derivative of the nucleon field to obtain field equations for if> and ip that are related by a 
simple adjungation. 

In constrast to the standard QHD Lagrangian density with Dirac 7 M matrices in the 
kinetic term, the 4x4 matrices 

T, = Y[(l + W)g, v + X, u ]-T, (5) 

are introduced in the DC model (g^ = diag(l, — 1, — 1, — 1) is the metric tensor). They 
contain quantities 

W = W + t-W T^ = % + r-% X ia , = X ia , + T-X ia , (6) 

which are parametrized up to quadratic terms in the isoscalar meson fields 

r (i) r (2) r (2o) 

= — o T m ° a + — r m » u 7 

m m 4 m 4 

r (i) r (») 

% = Xp V = -!±—Tr^jjifj3 v (8) 

m m 4 

and linear terms in the isovector meson fields 

r (i) r (i) 

W = ^6 % = ^p^ ^ = 0. (9) 
m m 

The fields W, T M , and X^ v describe the coupling of the meson fields to the derivative of 
the nucleon field. The term W is similar to the extension of the ZM model and leads to 



a density dependence of the effective self energies. As will be seen later, the terms 7^ and 
are the essential terms to give a momentum dependence to the scalar and vector self 
energies, respectively. Since the covariant derivative is used in (Jl|), there appear also higher 
order nonlinear couplings to the nucleon field in the Lagrangian density. 

In addition to the four coupling constants I\ in £ specifying the minimal nucleon-meson 
couplings, seven new couplings constants (Ta \ Icr , y2~\ T^ a \ T^ b \ F s , Tp ) appear. For 
the sake of convenience, mass factors are introduced to obtain pure numbers for the new 
coupling constants. The large number of new couplings makes the Lagrangian of the DC 
model very flexible. If the new coupling constants vanish the matrix fl5|) reduces to the 
usual 7 M matrix and we recover the original QHD Lagrangian density with only mininal 
couplings. In comparison to the widely employed models with non-linear selfcouplings of 
the cr-meson (a-meson and cj-meson) there are five (four) additional free parameters in the 
DC model. 

The non-linear or density- dependent RMF models discussed in the introduction are not a 
simple subcase of the present DC models. They could be included by making the couplings 
Ti density-dependent or replacing them by polynomials in a and uo^. We do not want ot 
make this further extension. As will be seen later, the higher order terms in eqs. (|7|) and 
produce similar effects as such additional terms. 

The field equations of nucleons and mesons are derived from the corresponding Euler- 
Lagrange equations. In the mean field approximation the meson fields are treated as classical 
fields by replacing them with their expectation values. Only positive energy states of the 
nucleons are taken into account in the calculation of the densities (no-sea approximation). 
Symmetries of the nuclear system can lead to a considerable simplification the equations in 
the actual calculation, e.g. in infinite nuclear matter. 

The field equation of the nucleons 

T^iD^ - (m - T a a - Y s f ■ tj $ + ~ (z<9%) ^ = (10) 
can be transformed to the usual Dirac form 

[ 7 ^-£ M )-(m-£)]V> = (11) 
by introducing the scalar self energy 

£ = S - T^id" (12) 
7 



and the vector self energy 



V» - Wid„ - X tlv id v 



with 



and 



s = r a a + r s r-5 + r M (rw + r p r ■ p) - -(z<9%) 



[(1 + W)g, u + + T p T-fT) - l -{id,W) - -(iPX^) 



(13) 



(14) 



(15) 



The self energies ([12]) and (|13|) in the DC model are differential operators that act on the 
nucleon field. This fact is the essential extension in the new model with derivative couplings. 
Thus the self energies describe an interaction that contains a state dependence in addition 
to the medium dependence from the mesons fields. Additionally, both S and depend on 
scalar and vector meson fields. 

From the field equation flTOp the continuity equation <9 M J M = for the current density 
operator = ipT^ip is immediately derived. Therefore, g = (J ) has to be interpreted as 
the conserved baryon density in the DC model. The usual vector density Qy = (^70^) is 
no longer a conserved quantity. Similar we find <9 M J M = for the isospin current density 
operator J M = tpT^Tip . 

The Lagrangian density £, leads to the field equations for the isoscalar mesons 

F (2) 

-L (7 



rrr 



1 + 2- 



r 



(2a) 



m 



r (26) 



4 c s 



rrr 



r (i) 

r CT & h ^ 

m 
r (i) 



(16) 
(17) 



and for the isovector mesons 



d^S + m 2 s 6 = T s g s + 



d^ + mijr = r p r + 



r (i) 

m 
r (i) 

m 



(18) 
(19) 



They contain additional terms as compared to the linear QHD model. The masses of the 
isoscalar mesons are no longer constant but become effectively density dependent. The 
minimal coupling of nucleons and mesons leads to source terms with the usual scalar densities 



Qs = (ipip) Qs = {ipTi/j) 



(20) 
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and the new current densities 

which depend on the usual vector current densities 

Jim = (V>7^) % = 



(21) 



(22) 



and the scalar densities (see appendix [A]). 

Additional source terms proportional to rS , r£ , Tf\ and appear in the field 
equations with the derivative scalar densities 

1 



D 



Q 



-D 



hi/ 



(23) 
(24) 



and the derivative current densities 





1 


1 D = 


2 




1 


f = 


2 



(25) 
(26) 



As shown, the scalar densities ( P3"| ) and (JPj) can also be expressed in terms of the derivative 
tensor densities 



t 



D 



L fJtV 



-( i/j'y^iDvTijj + (iD u ip)>y ll Tip ) 



(27) 
(28) 



which are not symmetric in the Lorentz indices. More explicit expressions for the densities 
can be found in appendix 



III. APPLICATION TO INFINITE NUCLEAR MATTER 

Infinite nuclear matter in its ground state is a homogeneous, isotropic and stationary 
system. Because of these symmetries, the set of field equations simplifies considerably. 
Densities do not depend on space-time coordinates and only the time-like component of the 
currents remains. Correspondingly, the meson fields are constant and the spatial components 
of the vector meson field vanish. Additionally, only the third component in isospin of the 
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isovector densities and isovector meson fields remains. The matrix is diagonal in isospin 
space and it is useful to introduce the abbreviations 

r ± = 7o(l + W ± + X 00± )-To ± T l± = ll (l + W±) 2 = 1,2,3 (29) 

for the time-like and space-like components of r M with the quantities 

W ± = W ± (VP) 3 T 0± = T ± (f ) 3 X 00± = X m ± (* 00 ) 3 ( 30) 

where the case of protons is denoted by + and neutrons by — , respectively. The momentum 
independent contributions to the self energies ([12]) and (|13|) are given by 

S ± = T a a ± T s (6) 3 + T 0± [r^ ± r p (p ) 3 ] (31) 

and 

V 0± = (1 + W± + X 00± ) [T^luo ± T p (p ) 3 ] . (32) 
The meson fields are directly obtained from the field equations ([TBI) - (|19D for given densities. 

A. Solutions of the generalized Dirac equation 

In order to calculate the various densities appearing in the field equations we first have 
to find the nucleon field from the the generalized Dirac equation with self-energies that are 
differential operators in the DC model. Solutions of the nucleon field equation are found 
with the plane wave ansatz 

^(p, a, r) = u(p, a, r) exp (-ip^) (33) 

for a nucleon with four momentum = (E,p) and positive energy E. Spin and isospin 
quantum numbers are specified by a and r. The scalar self energy 

S± = S± - T ±E (34) 

and the vector self energy 

Z ± = V 0± -(W ± + X 0Q± )E (35) 

for nucleons in nuclear matter are explicitly energy dependent in the DC model. They 
decrease with energy for positive T ± and (W± + Xqq±). Of course, this linear dependence 
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will give a reasonable description of realistic self energies only for not too high energies. The 
quantities T ±, W±, and X 0Q ± increase with increasing density and the energy dependence of 
the self energies becomes stronger. This behaviour also leads to a limitation in the density 
range where the DC model can be applied because of the increase of the Fermi momentum. 
Since the energy E already contains the rest mass m of the nucleon the self energies are 
rather insensitive to the nucleon energy if E < m. The energy dependence itself is density 
dependent and it is different for protons and neutrons if 1^ 7^ or 7^ 0. The ansatz 
fl3"3f) in the Dirac equation (|ll|) leads to the condition 

M T - E T a-p T 
a ■ p T -M T - E T 
for the spinor u(p, a, r) if the abbreviations 

M T = m-Sr + TorE (37) 

E T = (1 + W T + X 00r )E - V 0r (38) 
3 T = (1 + W T )P (39) 

are introduced. These quantities carry the index r = ± because the fields are not necessarily 
equal for protons and neutrons in asymmetric nuclear matter. The dispersion relation 

M 2 T =E 2 T -p 2 T (40) 



u(p, a, t) = (36) 



connects the energy 



and mass 



Er ^(AA + v^) (41) 



Mr = 



where the constants 



J^S (^r + B T V^+W) • (42) 

_ m- S T + B T Vor R _ T 0r . . 

^ ' ~VT = W~ " l + Wr + X^ (43) 

have been introduced. Not only the energy E T , but also the mass M T , depends on the 
momentum, the density and isospin. Obviously, reasonable solutions of the Dirac equation 
exist only for — 1 < B T < 1 and 1 + W T + X QOr ^ 0. The energy in the four-momentum of 
the nucleon in (RBI) is given by 



E T = (l + Wr+XoOr)- 1 



V 0t + 7 =L= (A t B t + v^f 



Pi 



(44) 
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Solutions of equation ( |36|) are found to be 



u(p,a,r) = J _ *° \C T (45) 

W 2[(1 + W T + X 00t )£ t - T 0r M r ] ^ £F^_ X(T 1 

with spin and isospin eigenfunctions Xa an d 6r- The normalization of these spinors is dis- 
cussed in appendix [B|. 

Negative energy solutions are obtained by a similar procedure.. Then the theory can be 
quantized in the usual way by expanding the nucleon field operator in terms of the complete 
set of solutions and imposing anticommutation relations for the creation and annihilation 
operators. However, in the mean field solution we do not need to carry out this standard 
procedure in detail here. 

B. Densities 

The densities that appear in the field equations are obtained by a summation over all 
occupied states in nuclear matter up to the Fermi momentum p^ for protons and neutrons, 
respectively. In the no-sea approximation only states of positive energy are considered. 
Proton and neutron constribution to the various densities can be obtained with the help of 
the projection operator P T = (1 + r(r) 3 )/2, e.g., the total baryon density 

Q = Jo = Q+ + Q- (46) 
and the third component of the isospin density 

'jo) =Q+~Q- (47) 



3 

are the sum and the difference of the proton and neutron densities 

Qt = §T Q P T il>) , (48) 

respectively. The other isoscalar and isovector densities are related to the corresponding 
proton and neutron densities in a similar manner. Considering the normalization of the so- 
lutions of the Dirac equation all densities can be expressed in terms of the three fundamental 
integrals 

~F 

II = [ T dp T p 2 T = \(p F T f (49) 

12 



1 2 



1 ?, 



Pi 



o dPr V^+M 2 



~F , rpF 

f r E F T -A\\^ ' L 



A r 



dPr pIVaJTW = ^fr x + ^-11 



(50) 
(51) 



with = (1 + W)pf and E F = y/A%~+ (Pr) 2 - The proton and neutron densities are then 
given by 

g T = (</Sr PW>> = 2^(1 + ^)~ 3 /[ = ^ (pf) 3 (52) 

with the spin degeneracy factor k = 2. This relation defines the corresponding Fermi 
momentum for a given nucleon density. Immediately, we obtain the scalar density 



QsT = (^p^) = c; 1 [B T r x + A T r 2 



(53) 



and the current density 



j 0r = (V7o^> = C~ l [I{ + A T B T F 2 



with 



C T = 7T 2 (1 + W T + X 00t )(1 + W T )\l - B 2 T ) 

With these results the relation 

g T = (1 + W T + X 00r )j 0T - Tq t q st 



(54) 



(55) 



(56) 



is easily confirmed. Considering the relation (|38| ) the derivative densities are most easily 
calculated from 



and 



where 



and 



ji = (4>E T P T ^) = (1 + W T + Xoo.)" 1 (gf T + V 0r g ST ) 



4 r = (^oErPrlp) = (1 + W T + Xo^y 1 ( Q ? + V 0t J 0t ) 



(57) 



(58) 



qI = (i>E T P T ^) = C-\l - B 2 )^ [A T (1 + Bl)r x + A\B T I T 2 + B T F 3 ] (59) 



q e t = (&y E+P T il>) = C;\l - Bl)^ [2A T B T I{ + A 2 T B 2 T F 2 + F 3 ] . 



(60) 
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The scalar derivative density is deduced from 

Qi = tU^" = (hoE T Pr*P) ~ (H " pP^) (61) 
= C - (1 + - <#) 

with (|B4j) , the dispersion relation (^0|), and 

e f T = §M T P T il>) = C~\\ - B 2 T y\ [2A T B T Il + A 2 T Il + B 2 I^] . (62) 

Nuclear matter is characterized by its equation of state, i.e. the energy per nucleon or 
pressure as a function of the nucleon density. The energy density is calculated from the 
energy-momentum tensor 



T x>1 = 



dC 



d{d x ip) 



d'V - g Xfl £ 



(63) 



where the sum runs over all fields (cp = ip,ilj,a,5,uj u , p v ). Since there are no contributions 
from the derivatives of the meson fields in nuclear matter the expectation value of the 
energy-momentum tensor reduces to 



(T x ^) = (v>rW) -g Xtx {£) 



with 



(64) 



(65) 



iir.^'^uj^ + m 2 p p^ ■ p p — m 2 cr 2 — m 2 5 ■ 5 . 

Only the meson fields contribute to {£) because the nucleon fields are solutions of the Dirac 
equation flTU|). With the help of equations (|3"7D, fl3"8|), and ( |B"4] ) the energy density 

(66) 



(Too) = ^T E ^) - <£> 

[lo(Er + V 0r ) - (M T - m + S T ) 

T 

[Qt - Qfr + (m- S T )g ST + V 0t j 0t ] 



and the pressure 



V 



2 2, 2-2 22 2r2 

m ^0 + m pP0 ~ m * a ~ m 8° 



\Y,(T tt ) = \^.p^) + (£) 



(67) 



i=i 



E, E M\ I 22, 2 S 22 2?2 
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are found. The DC model in the mean field approximation is thermodynamically consistent 
because the pressure ( |67| ) calculated from the energy-momentum tensor agrees with the 
thermodynamical pressure p t h = Q 2 -j- Q Furthermore the relation 



e + p = p + Q+ + P-Q- 



(68) 



is obtained with the chemical potentials 



p 7 



'1 + W T + X 00r ] 



V 0r + 



(a t b t + v^+WF 



(69) 



which are the energies 
Hove theorem [1551 holds. 



of the nucleons at the Fermi surface and the Hugenholtz-van 



Determing the properties of nuclear matter in the DC model for a given density g and 
neutron-proton asymmetry a = Qn ~ Qp is more involved than in standard RMF models. The 

Qn T" Qp 

meson fields have to be determined from the densities which themselves depend on the meson 
fields. A self-consistent solution is achieved by iteratively solving the set of equations until 
convergence is reached. 



C. Optical potential 

The Schrodinger equivalent optical potential V opt serves as a convenient means to charac- 
terize the in-medium interaction of a nucleon. It is obtained by recasting the Dirac equation 
(|TT|) into a Schrodinger-like equation for the large (upper) component of the nucleon spinor. 
In the DC model the momentum p in the generalized Dirac equation is multiplied by a factor 
(1 + W T ) in eq. (|39|) . In order to start for the non-relativistic reduction from the usual form 
of the time-independent Dirac equation 

[ 5 . p + p ( m _ E f )] 1> T ={E- ££) t/v (70) 

the effective self energies 

ff = S r + W T m = S T - T 0t E + W T m 

1 + W T 1 + W T [ ' 



and 



^cflf ^Or + W T E Vq t — X 00t E 



^" T+W T ~ 1 + W T {<Z) 
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are introduced. They depend linearly on the nucleon energy E if Tq t ^ and Xqq t 7^ 0. 
Since E contains the rest mass of the nucleon, there is only a small variation of the self- 
energies at small momenta of the nucleon. It is seen that the coupling vertices in the DC 
model with r£ , Tp, and lead to the energy-dependent scalar and vector self-energies, 
respectively. In the DC model the energy dependence results from the coupling of the 
Lorentz- vector meson fields to the derivative of the nucleon field. The Lorentz-scalar meson 
fields do not lead to an energy dependence of the effective self energies. The modification 
of the initial Lagrangian (HD by the term W alone does not introduce an energy dependence 
of the effective self energies. This resembles the ZM model in which the derivative coupling 
is a special case of the W contribution in the DC model. By a rescaling of the nucleon 
wave function in the ZM model, the derivative coupling was essentially removed, leading to 
a nonlinear coupling of the nucleon field to the cx-meson without an energy dependence of 
the self energies. 

The quantities T 0r and X QOt have different density dependencies. At low densities they 
are approximately proportional to the baryon density and the square of the baryon density, 
respectively. But with increasing density this proportionality is lost because of the nonlin- 
earity of the meson field equations. Even without an energy dependence the effective self 
energies ( [71] ) and (|7^) in the DC model show a nontrivial density dependence. 

There are various methods to derive the optical potential from the Dirac equation (170) ■ 
Here, the definition 



t t ycff ycff I ^ 



Sf) 2 -W) 2 (73) 



for the Schrodinger equivalent optical potential is used as in Refs. ||49|| . It has the property 
that V op tr rises linearly with the energy E of the nucleon if the self energies are energy 
independent. An alternative choice is the definition as used in Ref. |^5| which corresponds 
to multiplying (|73|) by m/E. In Ref. J5fj, |56| the optical potential is calculated from the 
difference of the energy E of the nucleon in the medium and the kinetic energy of a free nu- 
cleon without interaction at the same momentum p. This definition is more appropriate for 
a general relativistic definition of an optical potential but it is no longer a Schrodinger equiv- 
alent potential. It always leads to a constant optical potential for large nucleon momenta 
for energy independent self energies. For a meaningful comparison of optical potentials it 
is, of course, necessary to chose the same definition. 
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The optical potential (|73|) in the DC model is a quadratic function of the nucleon energy 
because the effective self energies ([H]) and (^) are linear functions of the energy. This 
behaviour is reasonable only in a limited range of energies and the model cannot be applied at 
arbitrarily high energies. But for nucleon energies up to about 1 GeV the linear dependence 
is a sufficient approximation. Indeed, self energies from Dirac phenomenology of proton 
scattering show an approximately linear decrease for small energies fJ5|. 

IV. PARAMETRIZATIONS AND PROPERTIES OF NUCLEAR MATTER 

The DC model introduces a large number of additional couplings as compared to the 
original RMF model with minimal meson-nucleon couplings and even with respect to usual 
non-linear models. The corresponding coupling constants have to be fixed by the limited 
number of nuclear matter properties. Here, we develop parametrizations within restricted 
versions of the DC model in order to investigate whether the model gives reasonable results. 

Symmetric nuclear matter is essentially characterized by the binding energy per nucleon 
E/A and the incompressibility K at saturation density g sat where the pressure p is zero. 
This leads to three conditions for the seven isoscalar coupling constants T a , T^, T w 
Tu\ T^ a \ and T^ b \ Furthermore, an adjustment of the relativistic effective mass M cS = 
m — S cfr at saturation density can be used to achieve a reasonable spin-orbit splitting in finite 
nuclei. A reduction of M eff increases the spin-orbit splitting. The four isovector coupling 
constants T$, Tg, T p} and in the DC model determine the symmetry energy J and its 
derivative with respect to the density which can be parametrized by L = 3g^. In all our 
parametrizations, the 5-meson is not considered, which is common practice in most RMF 
parametrizations, because its contribution seems not to be well determined in parameter fits 
to properties of finite nuclei. 

All our parametrizations of the DC models are fitted to give E/A = —16.0 MeV for the 
binding energy per nucleon at a saturation density of g sat = 0.150 fm~ 3 with an incom- 
pressibility of K = 240 MeV. These are reasonable values as compared to the properties 
of symmetric nuclear matter in other models. The parameters are adjusted to an effective 
nucleon mass at g sat of M cS = 0.6m close to the value of several relativistic mean field 
parametrizations with spin-orbit splittings of single-particle levels close to the experimental 
values. These four conditions require at least four isoscalar coupling constants in the model. 
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The assumemd symmetry energy of J = 32 MeV fixes one isovector coupling constant. 

In the first parametrization (DC1) only the usual minimal meson-nucleon couplings with 
r a , T w , and T p and the new couplings and are assumed to be non-zero. This leads 
to a model without energy-dependent effective self-energies because T M = and X^ v = 0. 
The non-zero quantity W contains only contributions quadratic in the cr-meson and w-meson 
fields and correspondingly, these mesons acquire a density-dependent effective mass. The 
DC1 parametrization has the same number of parameters as the standard non-linear models 
with self-interactions of the cr-meson. 

The second and third parameter set (DC2/DC3) are chosen to give a reasonable de- 
scription of the Schrodinger equivalent optical potential of the nucleon in symmetric nu- 
clear matter at saturation density g SSut . Since the optical potential extracted in Dirac phe- 
nomenology from proton-nucleus scattering becomes almost constant around the kinetic 
energy E — m = 1 GeV of the nucleon, we adjust the coupling constants to exhibit a maxi- 
mum of the potential at this energy. The DC model is not really applicable to higher energies 
because the empirical self-energies show a deviation from the linear energy-dependence and 
in our approach the optical potential is a quadratic function of the nucleon energy. The 
coupling constants r£r and have to be positive to generate a reduction of the self- 
energies with increasing nucleon energy. In order to keep the number of parameters small 
we require in this model the additional constraints + = and = ■ For 
small densities this condition cancels the density dependence of the u;-meson mass in eq. ( |T7D 
and makes the coupling coefficients to the derivative densities equal in eqs. (0,0). Now 
the five independent isoscalar parameters are uniquely determined by five conditions. 

In the parametrizations DC1 and DC2 the constant T p is fitted to give a symmetry energy 
of J = 32 MeV wheras Tp 1 ^ is assumed to be zero. The sets DC2 and DC3 share the same 
isoscalar coupling constants but in model DC3 the two coupling constants of the p-meson 
are adjusted to give both a symmetry energy of J = 32 MeV and a derivative L = 40 MeV 
of the symmetry energy at saturation density in symmetric nuclear matter. This value of 
L is typical for non-relativistic Skyrme Hartree-Fock models. An increase of leads to a 
decrease of T p in order to give the same symmetry energy. Precise values for the coupling 
constants in all DC models are presented in Table |. 

It is instructive to compare various properties of symmetric nuclear matter and neutron 
matter in the DC models with other typical relativistic mean field models, in particular with 
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TABLE I: Coupling constants in the DC models. The nucleon mass is m = 939 MeV and the 
meson masses are m a = 550 MeV, m u = 783 MeV, and m p = 763 MeV. 

model r g r p rg iff iff iff ij 2ffl) ^ 

DCl 10.996338 12.831929 3.783530 0.000 0.000 0.000 84.77 12.30 0.00 
DC2 10.536751 13.754212 3.917768 2.392 2.392 0.000 8.735 -211.54 211.54 
DC3 10.536751 13.754212 0.943337 2.392 2.392 5.525 8.735 -211.54 211.54 



standard non-linear parametrizations and models with density-dependent couplings. Here 



we chose the parametrization NL3 [29| which contains tertic and quartic self-interactions of 



the cr-meson field and the set TM1 [30] which takes also self-interactions of the u;-meson into 
account. The VDD model [13] assumes a dependence of meson- nucleon couplings on the 
vector density that generates so-called "rearrangement" contributions in the self-energies 
leading to a non-trivial density dependence. These parametrizations were obtained by fits 
to properties of finite nuclei and nuclear matter. 

In Figure p] the binding energy of symmetric nuclear matter is shown for different models 
as a function of the density. Below saturation density all parametrizations lead to almost 
identical results, but at higher densities there are characteristic differences. The DCl model 
without energy dependent effective self energies is very similar to the NL3 parametrization 
with the non-linear self-interaction of the a meson. The effective density-dependence of 
the meson masses in the DCl model with ri- 2 '' ^> T^ff apparently leads to an effective 
density dependence of the scalar and vector self-energies that is comparable to the standard 
non-linear model. This is easily observed by comparing the effective coupling constants in 
symmetric nuclear matter. They are extracted from the self-energies according to 

i i 

/ y \ 2 /v„ \ 5 



Tf=mA-J and T e « = ™ u ) . (74) 

The effective coupling constants as shown in Figures § and |3| are normalized to one at 
saturation density of the respective parametrization and they are independent of the meson 
masses. In the NL3 parametrization the effective coupling of the tu-meson is independent of 
the density and it is almost density independent in the DCl model. On the other hand, the 
effective coupling of the cr-meson shows in both models an increase for large densities. This 
behaviour is in contrast to all other RMF models. 
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FIG. 1: Binding energy per nucleon for symmetric nuclear matter as a function of the nucleon 
density in various relativistic mean field parametrizations. 

The TM1 parametrization achieves a slightly softer equation of state for symmetric nu- 
clear matter at high densities than the NL3 and DC1 models, respectively, and the effective to 
coupling decreases significantly with increasing density. In this model a quartic self-coupling 
of the cu-meson are taken explicitly into account in the Lagrangian density. It gives rise to 
a o>q contribution in the meson field equation that becomes important at high densities and 
consequently reduces the vector self-energy. 

The VDD model shows an even softer equation of state at high densities as compared 
to the above parametrizations. This behaviour is determined by the functional form of 
the density dependence of the coupling functions, which was chosen in order to describe 
effective coupling constants extracted from Dirac-Brueckner calculations of nuclear matter. 
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FIG. 2: Normalized effective coupling constant of the w-meson in symmetric nuclear matter as a 
function of the normalized nucleon density in various relativistic mean field parametrizations and 
Dirac-Brueckner calulations. 

In Figures |] und |3| the corresponding results of two DB calculations (BMB with Bonn 
B potential [fL8[| , tHM p2| ) are shown for comparison. Here, the strong decrease of the 



couplings with the density is evident although there is a considerable difference between the 
DB calculations. The DC2 and DC3 models, which are only different with respect to their 
isovector properties, also exhibit a reduction of the effective couplings similar to the VDD 
parametrization and the DB results. The derivative couplings in the Lagrangian generate 
an energy and density dependence of the self energies that leads to a significant decrease 
of the effective couplings. Consequently, the binding energy per nucleon rises much slower 
with density than in the NL3, DC1 and TM1 parametrizations. 
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FIG. 3: Normalized effective coupling constant of the cr-meson in symmetric nuclear matter as a 
function of the normalized nucleon density in various relativistic mean field parametrizations and 
Dirac-Brueckner calulations. 

The density dependence of the effective a-meson coupling and of the scalar self-energy is 
directly related to the decrease of the relativistic effective mass of the nucleon M eff which is 
shown in Figure ^ for symmetric nuclear matter. We recall that the effective mass of the DC 
models is fixed to a value of 0.6 at g sat . There a two distinctly different groups, especially 
at high densities. The parametrizations NL3 and DC1 show a very similar behaviour as 
for the equation of state. The effective mass of these sets is substantially smaller than the 
mass of the other parametrizations at high densities and the curve exhibits less curvature. 
In symmetric nuclear matter the effective nucleon masses in the DC2 and DC3 model are 
identical but depend except on the density of the medium also on the momentum of the 
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FIG. 4: Relativistic effective nucleoli mass in symmetric nuclear matter in units of the nucleon rest 
mass as a function of the nucleon density in various relativistic mean field parametrizations. 

nucleon. In Figure [| the effective mass of nucleons in the DC2/DC3 parametrizations is 
shown for the Fermi momentum p F . The effective mass in the VDD model is smaller than in 
the other models at saturation density which is reflected by the too large spin-orbit splitting 
in finite nuclei p3[ but the density dependence is more similar to the TMl, DC2 and DC3 



cases. 

In Fig. H we show the momentum dependence of the neutron effective mass of different 
models for different asymmetries. The effective nucleon mass M in the DC models is fixed 
to 0.6 m at saturation density g sat = 0.150 fm -3 for nucleons with the Fermi momentum 
in symmetric nuclear matter. The scalar self energy in the DC1 parametrization does not 
depend on the energy and therefore the effective mass is momentum independent. It increases 
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FIG. 5: Momentum dependence of the neutron effective mass in units of the nucleon mass for the 
DC parametrizations in symmetric nuclear matter (thin lines) and neutron matter (thick lines) at 
g = 0.150 fm -3 . (The curves for a = 0.0 in the DC2 and DC3 model are identical.) 



only slightly with the neutron-proton asymmetry a at constant density. In the DC2 and 
DC3 models M eS increases with the momentum and with the asymmetry. For a = the 
effective mass is identical for DC2 and DC3 but in asymmetric nuclear matter the absolute 
value and the momentum dependence of M eS are different. This effect is caused by the 
different values of T p and the non-zero r£° in the DC3 case. The momentum dependence 
becomes stonger at higher densities. 

The main difference between the DC2/DC3 parametrizations and the other models be- 
comes apparent when the Schrodinger equivalent optical potential is studied, which is shown 
in Fig. H for g = £> sat as a function of the nucleon energy. In all models without energy de- 
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FIG. 6: Dependence of the Schrodinger equivalent optical potential on the nucleon energy in 
symmetry nuclear matter at g = 0.150 fm~ 3 in various relativistic mean field models and in two 
fits (DPI, DP2) of proton- nucleus scattering in Dirac phenomenology. 

pendent self-energies the optical potential by construction rises linearly with the energy 
in contrast to the real part of the optical potential extracted from Dirac phenomenology 
f45| , [46]. This is clearly seen in Figure |[ At energies E — m below approx. 200 MeV the 
optical potential is attractive. It higher energies it becomes repulsive. The linear increase 
of the RMF models leads to a substantial overestimate of the optical potential at energies 
above a few hundred MeV. Only the DC2 and DC3 models show a significant reduction of 
the repulsive potential at high energies, even though it is still larger than the experimen- 
tal results from Dirac phenomenology. A fit of the parameters in the DC model with less 
constraints on the coupling constants than in the DC2/DC3 sets can certainly improve this 



25 




FIG. 7: Dependence of the Schrodinger equivalent optical potential on the nucleon momentum 
in symmetric nuclear matter at various densities for the parametrizations DC1 (thin lines) and 
DC2/DC3 (thick lines). The potentials at g = 0.5^ sa t (dotted lines), g = 1.0^ sat (solid lines), 
g = 1.5^ sat (dashed lines), and g = 2.0g sat (dot-dashed lines) are shifted downwards by 150 MeV, 
100 MeV, 50 MeV and MeV, respectively. Horizontal lines mark the range in the momentum 
where the potential crosses zero. 

agreement but one has to consider that the phenomenological scalar and vector potentials 
have an imaginary part which is absent in the DC model. Therefore, a comparison of the 
optical potentials at high energies has to viewed with caution. However, the DC model 
represents a significant qualitative improvement as compared to standard RMF models. 

The DC models also makes a prediction on the momentum dependence of the optical 
potential at different medium densities. In Figure [7|the results of the parametrization DC1 
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FIG. 8: Dependence of the Schrddinger equivalent optical potential on the density in symmetric 
nuclear matter for the parametrizations DC1 (thin lines) and DC2/DC3 (thick lines) at nucleon 
momenta of p = MeV (dotted lines), p = 200 MeV (solid lines), p = 400 MeV (dashed lines), 
and p = 600 MeV (dot-dashed lines). 

are compared to the DC2/DC3 parametrizations at four nuclear matter densities between one 
half and twice the saturation density. The DC1 parameter set exhibits the typical behaviour 
of a RMF model with energy-independent self-energies. The optical potential rises fast with 
increasing nucleon momentum and becomes repulsive at progressively smaller momenta as 
the density of the medium increases. The DC2/DC3 models with momentum-dependent self- 
energies shows a similar trend as the DC1 model at small densities but at higher densities 
the optical potential is much less repulsive than in the DC1 model. The shift in the zero of 
V opt is less pronounced and the nontrivial momentum dependence generates a considerable 
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FIG. 9: Dependence of the symmetry energy for symmetric nuclear matter on the nucleon density 
in various relativistic mean field parametrizations. 

curvature. The density-dependence of the optical potential for constant nucleon momentum 
is shown in Fig. |8] and exhibits significant differences when the set DC1 is compared to the 
sets DC2/DC3. The optical potential vanishes at zero density independent of the nucleon 
momentum. It becomes smaller with increasing density and reaches a minimum before it rises 
again. The minimum is deeper for smaller momenta. At lower densities the parametrizations 
DC1 and DC2/DC3 display a more or less similar behaviour but at high densities the optical 
potential in the model without momentum dependent self-energies strongly increases and 
becomes very repulsive whereas the slope of V opt in the other model is much smaller and the 
repulsion sets in much later. 

The various parametrizations also exhibit a clearly distinct behaviour in the density 
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FIG. 10: Binding energy per nucleon for neutron matter as a function of the neutron density in 
various relativistic mean field parametrizations. 

dependence of the symmetry energy which is presented in Figure |^. The sets NL3, TM1, 
DC1, DC2 show an almost linear increase of J over a large range of densities. In all these 
models the coupling of the nucleons to the p-meson is described by only one parameter, 
the constant V p of the minimal coupling. In the VDD model, the curve bends significantly 
with a smaller slope at saturation density because V p is a decreasing function with density. 
The effect is even larger in the DC3 parametrization where the symmetry energy derivative 
L at saturation density was fitted to a smaller but more realistic value than in the other 
relativistic models. Comparing only the parametrizations NL3, TM1, and VDD, which were 
fitted to properties of finite nuclei, one notices that the symmetry energy at saturation 
density of the NL3 and TM1 models is much larger than in the VDD model but that J is 
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very similar around a density of 0.11 fm~ 3 . 

The density dependence of the symmetry energy directly correlates with the slope of the 
equation of state for neutron matter which is shown in Figure |T0|. The parametrizations NL3, 
TM1, DC1, and DC2 show a much stiffer equation of state than the density-dependent VDD 
model and the DC3 model with the smaller symmetry energy derivative L. The difference 
of the models is especially apparent at small densities. In Refs. [55, ^] it was shown that 
the slope of the neutron matter equation of state is directly correlated with the neutron 
skin thickness of finite nuclei, i.e. the difference between neutron and proton radii. This 
relation holds in both nonrelativistic Skyrme Hartree-Fock models and relativistic mean- 
field models. Experiments to determine neutron radii of finite nuclei with higher precision 
than available today can lead to better constraints on the model parameters and to better 
predictions for the neutron matter equation of state. It can be expected that the DC3 model 
with the smallest slope will also give a small neutron skin thickness similar to several Skyrme 
parametrizations when applied to calculations of finite nuclei. 



V. SUMMARY AND OUTLOOK 



The DC model is an extension of standard relativistic quantum hadronic models in or- 
der to generate an effective momentum dependence of the self energies in the mean field 
approximation. This is achieved by introducing couplings of the meson fields to derivatives 
of the nucleon field in the Lagrangian density. The new contributions lead to additional 
source terms in the meson field equations and to a density dependence of the effective me- 
son masses. The effective interaction in the DC model is both medium and state dependent. 
The self energies in the Dirac equation are differential operators with a nontrivial density 
dependence. The effective mass depends in general on the density, the momentum and the 
neutron-proton asymmetry. 

We applied the DC model to nuclear matter and developed three parametrizations for the 
coupling constants under different conditions. The set DC1 without momentum dependent 
self energies is similar to previous nonlinear RMF models. It displays rather stiff equations 
of state for symmetric nuclear matter and neutron matter. The sets DC2 and DC3 lead 
to energy dependent self energies The equation of state becomes softer and is comparable 
to the RMF model with density dependent couplings. The Schrodinger equivalent optical 
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potential is much less repulsive than in standard RMF parametrizations at high nucleon 
energies. It shows an energy dependence more like optical potentials extracted from Dirac 
phenomenology. The DC model is also flexible to adjust isovector properties which affects 
the neutron matter equation of state. 

If the nucleon momenta are too large the assumption of a linear energy dependence of 
the self energies is no longer valid. In this case the DC model can no longer be applied. 
This deficiency could be improved by introducing couplings of the meson fields to higher 
derivatives of the nucleon fields. Since the medium density determines the Fermi momentum 
there is also a limitation of the model to not too large densities but the Fermi momentum 
increases only slowly with the density. 

Because of the limited number of constraints from nuclear matter properties the param- 
eter space of the DC model was not fully explored. It is worthwhile to apply the model 
to finite nuclei and to study the effects of the new coupling vertices in this case. Also 
important is the comparison to extensive experimental data on nucleon scattering. The 
parametrizations developed here can be used as a reasonable starting point for a further 
refinement. 

The approach of the DC model with fixed coupling constants but generalized interaction 
vertices in the Lagrangian density is in the spirit of non-linear RMF models. Alternatively, 
one can think of an extension of models with density dependent meson-nucleon vertices by 
introducing a dependence of the coupling functions on derivative densities. We currently 
also investigate the possibilities of such an approach and results will be reported in the 
future. 
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APPENDIX A: EXPLICIT EXPRESSIONS OF DENSITIES 

The conserved current densities in the field equations of the uj- and p-meson are given as 
J, = [(1 + W)g^ + X^f - % 6s + q s (Al) 
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(A2) 



and depend on the usual scalar densities and vector densities. The derivative scalar densities 



?? = A - - v„ • r = <r £ 



£ - r^f - W + iT p p, xf = gTt 



ID 

/IV 



(A3) 
(A4) 



with 



Q d s = 2(^7^- (id»i>h^}) = 9»Xv 

are obtained from a contraction of the derivative tensor densities 
1 ,r 



u 



t 



D 



2 
1 



^^D^ip + {iD v i\))-f^ )=t - TjJvjn - Y p p v ■ % 



(A5) 
(A6) 



(AT) 



( Tjr/niDvTi/; + (iD^^rip ) = & - Y^u v % - T p p u j p + iY p p u x j p (A8) 



where 



(A9) 
(A10) 



with the metric tensor g^ v . For the derivative current densities the expressions 



3^ = it - FuUnQs - T p p p ■ g s 

J,, = - VuU^Qs - Tpp^Qs + i^pp/j, x g s 



(All) 
(A12) 



with 



j [i 
3 d 



(A13) 
(A14) 



are found. 
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APPENDIX B: NORMALIZATION OF SPINORS 



The spinors fl45|) in the solutions of the Dirac equation are normalized according to 

u(p, a, t)T 0t u(p, a', r') = 5 aa >5 TT > (Bl) 

which guarantees that the zero component of the current density J M can be interpreted as 
the baryon density. Correspondingly, the relations 

u(p, cr, r)7 u(p, a', t') = T —^- — (B2) 

1 K ' (1 + W T + X 00T )E T - T 0T M T 



(1 + W T + X 00r )(l - Bl) \ ^fAj+M, 
and 

-t- \ (- i >\ M T 5 aa >5 TT > 

u{p,a,T)u(p,a ,t ) = ~ — (B3) 

(1 + W t + X 00t )E t -T 0t M t 

5vct'5 tt > B + ^ T 



(1 + W T + X 00r )(l -B*)\ VAJ+? T/ 

for the usual scalar and vector densities are found. In the calculation of the scalar derivative 
density also the relation 

~ 2 5 ,5 

u(p, a, tW • p T u(p, a', t') = r uu T l — (B4) 

(1 + W t + X 00t )E t -T 0t M t 

will be needed. 
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